Universal properties of infrared oscillator basis extrapolations 
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Recent work has shown that a finite harmonic oscillator basis in nuclear many-body calculations 
effectively imposes a hard-wall boundary condition in coordinate space, motivating infrared extrap- 
olation formulas for the energy and other observables. Here we further refine these formulas by 
studying two-body models and the deuteron. We accurately determine the box size as a function of 
the model space parameters, and compute scattering phase shifts in the harmonic oscillator basis. 
We show that the energy shift can be well approximated in terms of the asymptotic normaliza- 
tion coefficient and the bound-state momentum, discuss higher-order corrections for weakly bound 
systems, and illustrate this universal property using unitarily equivalent calculations of the deuteron. 
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I. INTRODUCTION 

Harmonic oscillator (HO) basis expansions are widely 
used in nuclear structure calculations, but limited com- 
putational resources often require that the basis be trun- 
cated before observables are fully converged. In such 
cases, a procedure to extrapolate results to infinite basis 
size is needed. Such schemes have conventionally been 
formulated using the basic parameters defining the os- 
cillator space, namely the maximum number of oscil- 
lator quanta TV and the frequency f2 of the oscillator 
wave functions. An alternative approach to extrapola- 
tions is motivated by effective field theory (EFT) and 
based instead on explicitly considering the infrared (IR) 
and ultraviolet (UV) cutoffs imposed by a finite oscil- 
lator basis [TJ |2J. This has recently led to a theoreti- 
cally motivated IR correction formula and an empirical 
UV correction formula [2j in which the basic extrapo- 
lation variables are an effective hard-wall radius L and 
the analogous cut-off in momentum, Ajjv- In terms of 
the oscillator length b = J h / (mfl) , rough estimates 
of these variables are L « y/2(N + 3/2)6 = L and 
A uv « V / 2(N + 3/2)fc/6 nog. 

The b dependence of L and Ajjv suggests that if the 
oscillator length is small enough (i.e., if the oscillator fre- 
quency is large enough), the UV correction will be neg- 
ligible compared to the IR correction. In this domain, 
an estimate for the energy in the truncated basis was de- 
rived in Ref. [2] based on an effective Dirichlet boundary 
condition at L: 
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where k^ = — 2mE co /h 2 is the binding momentum 
defined from the separation energy . Consideration of 
the tails of the HO wave functions motivated an improved 
choice for L given N and Ml [2J: 
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The extrapolation formula ([lj is the leading-order cor- 
rection to the ground-state energy once UV corrections 
can be neglected and once L exceeds the radius of the nu- 
cleus under consideration. Test calculations of few- and 
many-body nuclei using L = L' and with E^, A, and 
fcoo as fit parameters showed that the IR correction for- 
mula (fTl) can be used in practice [2]. (Note: The results 
in Refrpj were derived in the laboratory system with m 
the particle mass. Here for convenience we take m to be 
the reduced mass /z, which rescales b and fcoo but leaves 
the expressions unchanged.) 

In the present work we seek a more complete under- 
standing of this correction formula and to more accu- 
rately determine the hard- wall radius L. While the most 
useful application of Eq. (JlJ is to few- or many-body nu- 
clei, we specialize here to the two-particle case, which 
we can control and calculate precisely. In doing so we 
gain insight into the universal features of the IR extrap- 
olation, including its invariance to phase-shift equivalent 
potentials and its application to excited states. While the 
coefficient A was previously treated purely as a fit param- 
eter, we extend the derivation from Ref. [2] to show how 
it can be expressed in terms of the observables and 
the asymptotic normalization constant 7^, just as in the 
related Luscher-type formulas developed for lattice ap- 
plications [3 -6 \ . We examine the approximations leading 
to Eq. ([!]) and derive a corrected formula appropriate for 
weakly bound states, which is shown to work much better 
for the deuteron. 

Our strategy is to use a range of model potentials for 
which the Schrodinger equation can be solved analyti- 
cally or to any desired precision numerically to broadly 
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test and illustrate various features, and then turn to the 
deuteron with a set of phase-shift equivalent potentials 
for a real-world example. In particular we will consider: 



exact u(r) 
HO basis 



V sw {r) 
V cxp {r) 

V s (r) 
V q (r) 



-V Q 6{R-r) 



[square well] 
[exponential] 
[Gaussian] , 
[quartic] , 



(3) 
(4) 
(5) 
(6) 



where for each of the models we work in units with h = 1, 
reduced mass /u = 1, and express all lengths in units of 
R and all energies in units of h 2 //j,R 2 . For the realis- 
tic potential we use the Entem-Machleidt 500 MeV chi- 
ral EFT N 3 LO potential [7] and unitarily evolve it with 
the similarity renormalization group (SRG). These po- 
tentials provide a diverse set of tests for universal prop- 
erties. Because we can go to very high Ml and N for 
the two-particle bound states (and therefore large kuv)-, 
it is possible to always ensure that UV corrections are 
negligible. 

In Section [n] we determine a more accurate value for 
L than L' and show that the theoretically founded expo- 
nential form of the extrapolation is favored over Gaussian 
or power-law alternatives in practical applications. The 
accurate determination of the box radius L also allows us 
to compute scattering phase shifts directly in the oscil- 
lator basis. The derivation of the exponential form from 



Ref. [2] is extended in Section III to show that it depends 
only on observable quantities, and is therefore indepen- 
dent of the potential and has the same form for excited 
states. These formal conclusions are tested with model 
potentials and the deuteron with a realistic potential in 
Section [TV] In Section E] we summarize our conclusions 
and discuss the implications for applications to larger nu- 
clei. 



II. SPATIAL CUTOFF FROM HO BASIS 
TRUNCATIONS 

In this Section, we determine the spatial extent of a fi- 
nite HO basis. We start with empirical considerations be- 
fore presenting an analytical understanding. Finally, we 
use the knowledge of the spatial extent to compute phase 
shifts and demonstrate that the theoretically founded ex- 
ponential extrapolation law can be distinguished from 
other empirical choices. 



A. Empirical determination of L 

The derivation of the IR correction formula Eq. ([I]) 
in Ref. [5] starts from the observation that a truncated 
harmonic oscillator (HO) basis effectively acts at low en- 
ergies to impose a hard-wall boundary condition in coor- 
dinate space. In Fig. [T]we can see how this happens for a 
representative model case, a square well potential Eq. ^ 
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FIG. 1: (color online) (a) The exact radial wave function 
(dashed) for a square well Eq. (J3j with depth Vb = 4 (and 
H = fi = R=l)is compared to the wave function obtained 
from an HO basis truncated at N = 4 with hfl — 6 (solid). 
The spatial extent of the wave function obtained from the HO 
basis truncation is dictated by the square of HO wave function 
for the highest radial quantum number (dot-dashed), (b) The 
wave functions obtained from imposing a Dirichlet boundary 
condition at Lq, L' and L2 are compared to the wave function 
in truncated HO basis. 



with s-wave radial wave functions. In the top panel, the 
exact ground-state radial wave function (dashed) is com- 
pared to the solution in an oscillator basis truncated at 
N — 4 determined by diagonalization (solid). The trun- 
cated basis cuts off the tail of the exact wave function 
because the individual basis wave functions have a radial 
extent that depends on Ml (from the Gaussian part) and 
on the largest power of r (from the polynomial part). 
The latter is given by N = 2n + I. With N = 4 and 
I = 0, this means that n = 2 gives the largest power. 

The cutoff will then be determined by the n = 2 oscil- 
lator wave function, ?4i=2( r )j whose square (which is the 
relevant quantity) is also plotted in the top panel (dot- 
dashed). It is evident that the tail of the wave function 
in the truncated basis is fixed by this squared wave func- 
tion. The premise of Ref. [2] was that this cutoff is well 
modeled by a hard-wall (Dirichlet) boundary condition 
at r = L. If so, the question remains how best to quanti- 
tatively determine L given ./V and Ml. Before we present 
an analytical derivation of this quantity in the next Sub- 
section, we compare empirically L' from Eq. ^ and 



Li EE y/2(N + 3/2 + i)b 



(7) 



with integer i, which includes Lq as a special case. In the 
bottom panel of Fig. [T] we show the wave functions for 
several possible choices for L. L corresponds to choosing 
the classical turning point (i.e. the half-height point of 
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the tail of [u££ 2 (r)] 



'■); it is manifestly too small. Using 
Lq, which is the linear extrapolation from the slope at the 
half-height point, gives an improved estimate. However, 
choosing i = 2 (i.e., using L = L 2 = y/2(N + 3/2 + 2)6) 
is found to be the best choice in almost all examples. 




FIG. 2: (color online) Ground-state energies versus (a) Lo, (b) 
L'q, and (c) L2 for a Gaussian potential well Eq. Q with Vb = 
5 (and h = fj, = R = 1). The crosses are the energies from 
HO basis truncation. The energies obtained by numerically 
solving the Schrodinger equation with a Dirichlet boundary 
condition at L lie on the solid line. The horizontal dotted 
lines mark the exact energy, Eoo = —1.27. 

The most direct illustration of this conclusion comes 
from the bound-state energies. In the example in Fig. [T] 
the exact energy (in dimensionless units) is —1.51 while 
the result for the basis truncated at N = 4 is —1.33, 
which is therefore what we hope to reproduce. With Lq, 
the energy is —0.97, with L' Q it is —1.21, and with L 2 
it is —1.29. While this is only one example of a model 
problem, we have found that L 2 always gives a better 
energy estimate than L' (and L 3 is almost always worse) . 

Another signature that demonstrates the suitability of 
L 2 is that points from many different Ml and N values all 
lie on the same curve. Figures. [2] and [3] show the energies 
from a wide range of HO truncations for Lq, L'q and L 2 
for the Gaussian well and the square well potential, re- 
spectively. The energies for different Ml and N lie on the 
same smooth and unbroken curve if we use L 2 but not 
with the other choices. For L = Lq and L = L' , one finds 
that sets of points with different Ml but same N fall on 
smooth, iV-dependent curves. For the square well, there 
are small discontinuities visible even for L = L 2 . At the 
square well radius, the wave function's second derivative 
is not smooth, and this is difficult to approximate with a 
finite set of oscillator functions. This lack of UV conver- 
gence is likely the origin of the very small discontinuities. 
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FIG. 3: (color online) Ground-state energies versus (a) Lq, 
(b) L' , and (c) L2 for a square well potential well Eq. Q 
with Vb = 4 (and h = fj, = R = 1). The crosses are the 
energies from HO basis truncation. The energies obtained by 
numerically solving the Schrodinger equation with a Dirichlet 
boundary condition at L lie on the solid line. The horizontal 
dotted lines mark the exact energy, Eoo = —1.51. 



As a further test, we solve the Schrodinger equation with 
a vanishing Dirichlet boundary condition (solid lines in 
Figs. [2] and [3]) and compare to the energies obtained from 
the HO truncations (crosses). The finite oscillator basis 
energies are well approximated by a Dirichlet boundary 
condition with a mapping from the oscillator Ml and N 
to an equivalent length given by L 2 . Note that for large 
N , the differences between Lq, L'q and L 2 may be smaller 
than other uncertainties involved in nuclear calculations, 
but for practical calculations one will want to use small 
N results, where these considerations are very relevant. 



These results from model calculations are consis- 
tent with those from realistic potentials applied to the 
deuteron. To illustrate this, we use the N 3 LO 500 MeV 
potential of Entem and Machleidt [7]. We generate re- 
sults in an HO basis with Ml ranging from 1 to 100 MeV 
and N from 4 to 100 (in steps of 4 to avoid HO artifacts 
for the deuteron [5j). We then restrict the data to where 
UV corrections are negligible (see Section IV C ). Figure [4] 
shows that the criterion of a continuous curve with the 
smallest spread of points clearly favors L 2 . Similar com- 
ments apply to the computation of the radius. Figure [5] 
shows that the numerical results for the squared radius, 
when plotted as a function of L 2 (but not as a function 
of Lq), fall on a continuous curve with minimal spread. 
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FIG. 4: (color online) Ground-state energies versus (a) Lo, 
(b) L' , and (c) i 2 for the Entem-Machleidt 500 MeV N 3 LO 
potential [7]. The horizontal dotted lines mark the exact en- 
ergy, Eev = -2.2246 MeV. 
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FIG. 5: (color online) Deuteron radius squared versus (a) Lq 
and (b) L 2 for the Entem-Machleidt 500 MeV N 3 LO poten- 
tial [7]. The horizontal dotted lines mark the exact radius 
squared, r 2 ^ = 3.9006 fm 2 . The insets show a magnification 
of data at smaller lengths L n . 



B. Analytical derivation of Li 

Naturally, the squared momentum operator p 2 is the 
key for understanding the IR properties of the harmonic 
oscillator basis. Let us start with the spectrum of p 2 in 
the oscillator basis. In a finite basis with energies up to 
E = (N + 3/2)M7, the operator p 2 must be viewed as 
p 2 Q(E — p 2 /(2m) — (m/2)Q 2 r 2 ), where denotes the 



unit step function. Let us compute the number M(k) of 
s-wave states up to a momentum k as a first step. We 
find 



M(k) = Tr 



ir)(h 2 k 2 -p 2 )Q[E-^-- r ^tfr 2 



hk 



JL J dp J dr e(n 2 k 2 -p 2 



-hk 



x 6 [E 



2m 



-Q 2 r 2 



(8) 



the semiclassical approximation and 
as a phase-space integral. We as- 



Here, we apply 
write the trace 
sume h 2 k 2 / (2m) < E, perform the integrations and use 
E/(Ml) = N + 3/2. This yields 



M(k) 



bk 
2^ 



y/2N + 3-b 2 k 2 

TV + 3/2 . bk 
arcsm 



V2N + 3 



(9) 



where b is the oscillator length. Figure [6] shows a compar- 
ison between the quantum mechanical staircase function 
and the semiclassical estimatej9| for N = 32. For suffi- 
ciently small values of kb <C V 2N, the number of s-wave 
momentum eigenstates grows linearly, and inspection of 
Eq. ([9]) shows that the slope at the origin is Lo/n semi- 
classically. The linear growth of the number of eigen- 
states of p 2 with k clearly demonstrate that — at not 
too large values of kb — the spectrum of p 2 in the oscil- 
lator basis is indistinguishable from the spectrum of p 2 in 
a spherical box. For the determination of the box radius 
L, we note that the lowest eigenvalue of p 2 is (ir/L) 2 . 




FIG. 6: (color online) The staircase function of the s states 
of the operator p 2 in a finite oscillator basis with N = 32 
(black) compared to its semiclassical estimate (smooth red 
curve). M(k) denotes the number of states of the operator p 2 
with eigenvalues p 2 < h 2 k 2 . 

In what follows, we analytically compute the smallest 
eigenvalue K^ in of p 2 in a finite oscillator basis and will 
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see that K m i„ = ir/L 2 . In the remainder of this Sub- 
section, we set the oscillator length to one. We focus on 
s-waves and thus consider wave functions that are regular 
at the origin, i.e. the radial wave functions are identical 
to the odd wave functions of the one-dimensional har- 
monic oscillator. 

The localized eigenfunction of the operator p 2 with 
smallest eigenvalue k 2 is 



sin nr 
, 



0<r<l 



We employ the s-wave oscillator functions 



(10) 



¥> 2 „+i(r) = (-in 



2nl 



T(n + 3/2) 

-1/2 



LI (r 2 ) 



2\ Sr 



(7rh 2,l (2n + l)!) 



H 



2n+l 



(r)e 



1/2 

with energy E = (2n + 3/2)hQ. Here, L n ' denotes the 
Laguerre polynomial, and it is convenient to rewrite this 
function in terms of the Hermite polynomial H n . We 
expand the eigenfunction ( 10 ) as 



^2 ° 2n ^ 

n=0 



(11) 



Before we turn to the computation of the expansion coef- 
ficients C2ti+i(k), we consider the eigenvalue problem for 
the operator p 2 . We have 



p 2 = a) a 



( a ty 



(12) 



where a and a) denote the annihilation and creation op- 
erator for the one-dimensional harmonic oscillator, re- 
spectively. The matrix of p 2 is tridiagonal in the os- 
cillator basis. For the matrix representation, we order 
the basis states as ((pi, tp3, ip$, . . .). Thus, the eigenvalue 



problem p 



becomes a set of rows of coupled 



linear equations. In an infinite basis, the eigenvector 
(ci(re), C3(/t), cg(K), . . .) identically satisfies every row of 
the eigenvalue problem for any value of k. In a finite 
basis ((pi, (f3, cps, . . . ip 2n +\), however, the last row of the 
eigenvalue problem 



(2n + 3/2-K 2 )c 2r , 



1 



-i( K ) = -V2nV2n + 1c 2 „_i(k) 



,( 13 ) 

can only be fulfilled for certain values of k, and this is the 
quantization condition. To solve this eigenvalue prob- 
lem we need expressions for the expansion coefficients 
C2ri+i(^) for n>l. Those can be derived analytically as 
follows. 



We rewrite the eigenfunction ( 10 1 as a Fourier trans- 
form 



dkip K (k) sin kr 



(14) 



and expand the sine function in terms of oscillator func- 
tions as 

I — oo 

sinfcr = , r- ^2(-l) n (p 2n+ i(r)<p 2n+1 (k) . (15) 



Thus, the expansion coefficients in Eq. (11) are given in 
terms of the Fourier transform tfi K (k) as 



C2n+i(«) = (-!)"/ dkr)) K (k) V i n+ x{k) . (16) 



So far, all manipulations have been exact. We need an 
expression for c 2n +i(K) for n» 1 and use the asymptotic 
expansion 

^ n+ i (k) « ^r^' sin( Vi^+Sfc) , (17) 



which is valid for |fc| <C V2n, see [9]. Using this approx- 
imation, one finds (making use of Fourier transforms) 



VT2n)! 
, /4 (2n-l)H . 



sin(V4n + 3k) , (18) 



with k < ?r/V4n + 3 due to Eq. fll0| . 

Let us return to the solution of the quantization con- 
dition (13 1. We make the ansatz 



~2A ' 



(19) 



and must assume that A > 0. This ansatz is well mo- 
tivated, since the naive semiclassical estimate k — tt/Lq 
yields A = 0. We insert the expansion coefficients ( 18 ) 



into the quantization condition ( 13 ) and consider its 



leading-order approximation for n>l and n ^> A. This 
yields 



A = 2 



(20) 



as the solution. Recalling that a truncation of the basis 
at y>2n+i corresponds to the maximum energy E = (2n + 
3/2)Ml, we see that we must identify N = 2n. Thus, 
«min = k / L 2 is the lowest momentum (or minimum step 
of momentum) in a finite oscillator basis with n 1 
basis states (and not l/b as stated in Ref. [I]). It is 
clear from its very definition that tt/L 2 is also (a very 
precise approximation of) the natural infrared cutoff in a 
finite oscillator basis, and that L 2 (and not b as stated in 
Refs. PZHH3]) is the radial extent of the oscillator basis 
and the analogue to the extent of the lattice in lattice 
computations j3]. 

The derivation of our key result /t m i n = tt/L 2 is based 
on the assumption that the number of shells N fulfills 
N > 1. Table [I] shows a comparison of numerical results 
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0.5310 


0.5642 


16 


0.5035 


0.5031 
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0.4585 


0.4582 


0.4791 



TABLE I: Comparison between the lowest momentum Kmin, 
7r/Z/2, and tt/Lq for model spaces with up to N oscillator 
quanta. 



for /c min in different model spaces. We see that ir/L 2 is 
a very good approximation already for N = 2, with a 
deviation of about 1%. 

Note that this approach can be generalized to other lo- 
calized bases. As the number of basis states is increased, 
the (numerical) computation of the lowest eigenvalue of 
the momentum operator p 2 yields the box size L corre- 
sponding to the employed Hilbert space, and results can 
then be extrapolated according to Eq. ([I]). 



C. Scattering phase shifts 

The argument for computing scattering phase shifts is 
as follows: The oscillator basis appears as a spherical 
box of size L. For low momenta we have L = L 2 , but 
at higher momentum L deviates slightly from L 2 , and 
can be determined from the eigenvalues of the operator 
p 2 . Thus, the positive-energy states computed in the 
oscillator basis can be used to extract phase shifts. 

In a fixed harmonic oscillator basis (N,MY), the com- 
putation of the phase shifts for a given partial wave 
2S+1 lj with orbital angular momentum I proceeds as fol- 
lows: First, one computes the discrete eigenvalues p 2 of 
the operator p 2 for orbital angular momentum I. Second, 
we need to determine the momentum dependent box size 



Assuming that the 

■th 



momentum eigen- 

genstate of a spher: 
zero of the spherical Bessel function. 
determines L(pi). We evaluate the 



Li = L(p, 

state is the i th eigenstate of a spherical box, we must 
determine the i 
Thus jiipiLi/H) 
smooth function L(p) for arbitrary momentum p by in- 
terpolating between the discrete momenta p t . Third, we 
compute the discrete positive energies = fi 2 fc|/(2m) 
of the neutron-proton system in relative coordinates for 
the partial wave 2S+1 lj, and compute the phase shifts 
from the Dirichlet boundary condition at r = L, i.e. 



Here iji is the spherical Neumann function. In practice 
one repeats this procedure for several values of Ml in 
order to get sufficiently many datapoints that fall onto a 
smooth curve. 

As examples we compute the scattering phase shifts 
for the 1 So and 3 Pi partial waves in model spaces with 
N = 32 and Ml = 20, 22, . . . , 40 MeV. Our calculations 
are based on the Entem-Machleidt 500 MeV chiral EFT 
N 3 LO potential [7J. Figures [7] and [8] show the results and 
compares them to the numerically exact phase shifts. For 
smaller N than our current choice, the computed phase 
shifts start to deviate from exact phase shifts at higher 
energies. However, if one is interested only in low-energy 
phase shifts and observables such as the scattering length 
and the effective range, a smaller harmonic oscillator ba- 
sis is sufficient. 
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FIG. 7: (color online) The 1 So phase shifts (in degrees) of the 
N 3 LO chiral interaction (solid line) compared to the phase 
shifts computed directly in the harmonic oscillator basis (cir- 
cles). 
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FIG. 8: (color online) The 3 Pi phase shifts (in degrees) of the 
N 3 LO chiral interaction (solid line) compared to the phase 
shifts computed directly in the harmonic oscillator basis (cir- 
cles). 
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There are other methods to compute scattering phase 
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shifts in the harmonic oscillator basis. Bang et al. [TT] 
used the method of harmonic oscillator representation of 
scattering equations (HORSE) for this purpose, and more 
recent works [T3] H3] computed phase shifts to develop an 
EFT for nuclear interactions directly in the oscillator ba- 
sis [10] . References [12l [13] build on the results by Busch 
et al. [Q] and their generalization [TS] to finite range 
corrections, and extract scattering information from the 
energy shifts of bound states in a harmonic oscillator po- 
tential. The resulting EFTs are quite efficient for con- 
tact interactions and systems such as ultracold trapped 
fermions, but nuclear potentials with a finite range re- 
quire an extrapolation of fl — > [T3] . The approach pre- 
sented in this Subsection is more direct, as no external 
oscillator potential is employed. We note that the anal- 
ysis presented in this Subsection can easily be extended 
to coupled channels as well. 

Finally, we note again that the approach of this Sub- 
section can be utilized in other localized basis sets. All 
that is required is the diagonalization of the operator p 2 
in the employed basis set, which yields the (momentum 
dependent) box size. 



D. Functional dependence of extrapolation 

The extrapolation formula ([I]) with L = L2 is theo- 
retically founded. How well can the specific form of this 
extrapolation be distinguished from other popular empir- 
ical choices? To address this question, we test possible 
functional dependences of the energy correction AE on 
L. The most common extrapolation schemes employ an 
exponential in N (or equivalently a Gaussian dependence 
on L), 



E(N) = E 



C 



N e 



-b N N 



(22) 



where Cn and arc determined separately for each Ml 
(with the option of a constrained fit of a common Eoo for 
special Ml values). Thus, unlike the extrapolation based 
on L, there is no universal variable and no distinction 
between IR and UV regions in Ml. However, empirically 
the form in Eq. ( 22 ) seems to work quite well [5J [ToTJT§] . 
Recently, Tolle et al. [2U] investigated the convergence 
properties of genuine and smeared contact interactions in 
an effective theory of trapped bosons and found that the 
smearing changed a power law dependence of the conver- 
gence to an exponential dependence. Here we will con- 
sider all three functional dependences on L: exponential, 
Gaussian, and power law. 

A purely empirical test can be made for our models and 
the deuteron because we can calculate the exact Eoo , plot 
&E(Li2) = EiLi) — Eoo against L2, and then attempt to 
fit each of the three choices of AE(L 2 ). Figure [9] shows 
the results for a representative model potential (a Gaus- 
sian) with moderate depth while Fig. 10 shows the results 
for the deuteron. The plots are made so that the candi- 
date form would yield a straight line if followed precisely. 
We see that the exponential form is an excellent fit for 




FIG. 9: (color online) The IR energy correction AEl versus 
Z/2 for a Gaussian potential well Eq. |5| with Vo = 5 (and 
fi — fi = R = 1) using a wide range of N and HQ,. The 
energies are fitted with (a) exponential, (b) Gaussian, and (c) 
power law dependence on L^. 



the model throughout the range of L2 and a reasonable 
but not perfect fit for the deuteron. For the deuteron, the 
weak binding is a challenge as it requires very large values 
of L2 for extrapolations. Corrections to weak binding will 
be derived in Section |ITt[ In contrast to the exponential 
extrapolation, Gaussian and power law fits fail over the 
full range of L 2 . This is consistent with Tolle et al. [20j . 
For limited ranges of L 2 a Gaussian does provide a rea- 
sonable fit (and should give a good extrapolation for Eoo 
if close enough to convergence) , but not globally. 

At this stage we have empirically verified the useful- 
ness of the extrapolation ([I]) in a very controlled setting. 
This corroborates the study in Ref. [2] and applications 
in Rcfs. [3JJE2]. The fit result for koo has generally been 
quantitatively consistent with nucleon separation ener- 
gies (note, however, the case of 6 He in Ref. [2]), but the 
constant A was not identified with physical quantities. 
The next section will express A in terms of observables 
for the two-particle system and present corrections to the 
extrapolation law ([lj. 



III. UNIVERSAL FORMULAS FOR IR 
CORRECTIONS 

In this Section we revisit the derivation of Eq. and 
obtain an expression for the coefficient A in terms of the 
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FIG. 10: (color online) The IR energy correction AEl versus 
1/2 for the deuteron calculated with the chiral EFT potential 
from Ref. [7] using a wide range of TV and HQ. The energies 
are fitted with (a) exponential, (b) Gaussian, and (c) power 
law dependence on L-z- 



bound-state asymptotic normalization coefficient (ANC) 
7oo and . This is in close analogy to correction formu- 
las for energies calculated with lattice regularization for 
periodic and hard wall boundary conditions j3H6j. Be- 
cause fcoo and are measurable, the result is universal 
in the sense that it is the same for any potential that 
reproduces the experimental observables for the bound 
state. The parameters in Eq. ([I]) can be fully predicted 
and tested against precise numerical fits for both our 
models and the deuteron, which is carried out in Sec- 
IV Corrections to Eq. ([I]) derived below are found 



tion 

to be quantitatively important for shallow bound states. 



A. Linear energy approximation 

Our first approximation to the IR correction is based 
on what is known in quantum chemistry as the linear en- 
ergy method |23j . Given a hard- wall boundary condition 
at r = L beyond the range of the potential, we write the 
energy compared to that for L = oo as 



Er 



AE r 



(23) 



We seek an estimate for AEl, which is assumed to be 
small, based on an expansion of the wave function in 
AE^. Let %(r) be a radial solution with regular bound- 
ary condition at the origin and energy E. For convenience 



in using standard quantum scattering formalism below, 
we choose the normalization corresponding to what is 
called the "regular solution" in Ref. [24] . which means 
that ue(0) = and the slope at the origin is unity for all 
E. We denote the particular solutions UE L (r) = Uz,{r) 
and ue x {t) = Uoo(r). Then there is a smooth expansion 
of ue about E = E^ at fixed r, so we approximate |23j 



u L (r) sa Uoo(f) + AE L 



du E {r) 



dE 



O(AEl) , 



(24) 



for r < L. By evaluating Eq. (24) at r — L with the 



boundary condition ul{L) = 0, wc find 

du E (L) 



AEr 



-Uoo(i) 



dE 



(25) 



which is the estimate for the IR correction. 

We can check the accuracy of the linear energy ap- 
proximation ( 24 ) by numerically solving the Schrodinger 



equation with a specified energy. This determines L as 
the radius at which the resulting wave function vanishes. 
Then we compare this wave function for r < L to the 
right side of Eq. (|24[ ) , with the derivative calculated nu- 
merically. Figure [TTj shows representative examples for a 
deep and shallow Gaussian potential. In these examples 
and other cases, the approximation to the wave function 
is good, particularly in the interior. The estimates for 
AE L using the right side of Eq. ( [25] ) are within a few to 
ten percent: 0.68 versus 0.70 and 0.050 versus 0.055 for 
the two cases. 
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FIG. 11: (co lor online) Testing the linear energy approxima- 
tion Eq. (24) for (a) deep (Vo = 10) and (b) shallow (Vo = 2) 
Gaussian potential well Eq. (5| (h = fj, = R = 1). The 
solid lines are the exact solutions Uc(r) for energies —3.5 and 
—0.020, respectively, whose zero crossings determine the cor- 
responding values for L. 
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The good approximation to the wave function suggests 
that for the calculation of other observables the linear en- 
ergy approximation will be useful. For observables most 
sensitive to the long distance (outer) part of the wave 
function, such as the radius, this has already been shown 
to be true [2]. But the good approximation to the wave 
function at small r means that corrections for short-range 
observables should also be controlled, with the dominant 
contribution in an extrapolation formula coming from the 
normalization. 



B. Complete IR scaling 

Next we derive an expression for the derivative in 
Eq. ( p5j ) . We assume we have a single partial- wave chan- 
nel and reserve the generalization to coupled channels 
(e.g., for a complete treatment of the deuteron) for fu- 
ture work. For general E < 0, the asymptotic form of 
the radial wave function for r greater than the range of 
the potential is (using the notation of Ref. [2]) 



with 



u E [r) — > A E {e 
,(r) A^e-^ 



k E r 



+k E r 



(26) 



We take the 



for E = Eoo. 

derivative of Eq. ( 26 1 with respect to energy, evaluate at 
E = using oiE^a = an d dkE/dE = —fi/(Jl < &e), to 
find 



du E (r) 
dE 



A 



daE 
° ~dE 

dA E 



+ fco 



A r 



/i r 



dE 



We now evaluate at r 
term dominates: 



L and anticipate that the e 



(27) 



+k x L 



du E (L) 



dE 



A* 



dar 



dE 



C>(e 



-fcooM 



Substituting Eq. (|28|) into Eq. (25 1, we obtain 

AE, = - — ± e- 2k °° L + 0{<r Ak ~ L ) 



daE 




dE 





(28) 



(29) 



which is in the form of Eq. ([!]). Note that this result is 
independent of the normalization of the wave function. 

To calculate the derivative explicitly, we turn to scat- 
tering theory, following the notation and discussion in 
Ref. |24) . In particular, the asymptotic form of the regu- 
lar scattering wave function cf>i_ k for orbital angular mo- 
mentum I and for positive energy E = h 2 k 2 /2/j, is given 
in terms of the Jost function fz(fc) |24j . 



Mk)K(kr)-U(-k)h+(kr)} 



(30) 



where the hf functions (related to Hankel functions) be- 
have asymptotically as 



hf (kr) 



t2= e ±i{kr-l7r/2) 



(31) 



The ratio of the Jost functions appearing in Eq. ( 30 1 gives 
the partial wave S'-matrix Si(k): 



si(k) = 



fi(-k) 



(32) 



which is in turn related to the partial-wave scattering 
amplitude fi(k) by 



fi(k) 



si(k) - 1 
2ik 



(33) 



We will restrict ourselves to I = for simplicity; the 
generalization to higher I is straightforward. 

To apply Eq. (30) to negative energies, we analytically 



continue from real to (positive) imaginary k. So, 



(t>o,ik E ( r ) 



{U^ E )e kEr - f (-ik E )e- kEr ) 



, — k E r 



hi-iks) c k Er \ 
fo(ik E ) ' 



(34) 



where R is the range of the potential. Upon comparing 
to Eq. ( 26 ) we conclude that 



as 



fo(ik E ) 
fo(-ik E ) 



1 



s (ik E ) 



(35) 



Note that Eq. ( 35 1 is consistent with the bound-state 
limit of Eq. ( 26 1 : at a bound state where E^ = 



there is a simple pole in the S matrix, which 
means ue = as expected (no exponentially rising 
piece). 

From Ref. [21] we learn that the residue as a function 
of E of the partial wave amplitude fi(E) at the bound- 
state pole is (-l) i+1 7^ft 2 /2^, where is the ANC. The 
ANC is defined by the large-r behavior of the normalized 
bound-state wave function: 

Wnor m (r) r -^ 7 ooe-^ r . (36) 

Thus, near the bound-state pole (with E — h 2 k 2 /2/i), 



fo(k) 



/ OO 



2 f i(E-E 00 ) k 2 + k^ 



or, using Eqs. (33 1 and (35), 



Now, 



a E {k) 



daE 
ItE 



kl 



k 2 + k lo ~ 2i &7, 



da E /dk\ k=iko 
dE I dk\ k=ikoc 



so using Eq. ( 38 1 we find 

daE 



dk 



(37) 



(38) 



(39) 



(40) 
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and therefore 



IV. TESTS OF IR CORRECTION FORMULAS 



da i 



dE 



(41) 



Putting it all together, we have 

II «oo7oo „-2fc oc ,I 



AEr 



/' 



0(e 



) , (42) 



in agreement with Eq. (JlJ, but now we have identified 
A = rPkoai^/n. 

If we apply this correction for a weakly bound state, 
such that koo is small, we may not be justified in neglect- 
ing the second term on the right side of Eq. (27). If we 
keep it instead, then Eq. (28) becomes 



du E (L) 
dE 



A p +k x L 




rah 



and we have a modified infrared scaling given by 



(AE L ) r 



M (1 - lloLe 



(43) 



(44) 



We will test both Eqs. d42h and (44) in Section IV 



C. Relation to Liischer-type formulas 



In this Section we test direct fits of Eq. ([IJ, which 
has three parameters, and the specialized expressions for 
AEl in Eqs. (42 ) and (44), which have no free parameters 
if we take fcoo and 7oo from the exact solutions. Based on 
the results presented in Sect, [ilj we use L2 in all our fur- 
ther analyses. It is important that we isolate the IR cor- 
rections in making these tests. The truncation in the HO 
basis also introduces an ultraviolet error inversely propor- 
tional to the ultraviolet cutoff Auv ~ y // 2fihfl(N + 3/2). 
In the results here we use combinations of Ml and N val- 
ues such that the UV error in each case can be neglected 
compared to the IR error. (This is verified quantitatively 
by using a fit ansatz from Ref. [2] for the UV correction, 
which is assumed to be independent of the IR correction.) 

For each of the model potentials, the radial 
Schrodinger equation is accurately solved numerically in 
coordinate space for the energy, which yields k^, and 
the wave functions. The asymptotic normalization coef- 
ficient 700 is found by multiplying the wave function by 
gfcoor anc j reac | m g ff it s asymptotic value. This is illus- 
trated in the inset of Fig. [l2j which also shows the onset 
of the plateau that defines the asymptotic region in Li 
where we expect our correction formulas to hold. For the 
deuteron, the Hamiltonian is diagonalized in momentum 
space to find fcoo, and then an extrapolation to the pole 
is used to find the s-wave and d-wave ANCs [25]. In 
the present work we use only the s-wave ANC for the 
deuteron. 



Starting with the seminal work of Liischer [5] , a wide 
variety of formulas have been derived for the energy 
shift of bound states in finite- volume lattice calculations. 
The usual application is to simulations that use periodic 
boundary conditions in cubic boxes (e.g., see Ref. [5]). 
The recent work by Pine and Lee 01 [5] extend the deriva- 
tion to hard- wall boundary conditions using effective field 
theory for zero-range interactions and the method of im- 
ages. The result for AEl in a three-dimensional cubic 
box has a different functional form than found here (the 
leading exponential is multiplied by 1 / L with that geom- 
etry) and the subleading corrections are parametrically 
larger. 

However, because the HO truncation we consider is 
in partial waves, the one-dimensional analysis and for- 
mula from Ref. 5\ are applicable (because fcoo and 700 are 
asymptotic quantities, the result for zero-range interac- 
tion is actually general for short-range interactions). The 
method of images can be applied in a one-dimensional 
box of size 2L after specializing to a particular partial 
wave and then extending the space to odd solutions in r 
from —00 to +00. The leading-order finite- volume cor- 
rection agrees with Eq. (42), and the first omitted term 
is of the same order. The methods presented in [11 [5] can 



A. Universal properties 



be used to extend the present formulas to higher orders 
and more general cases, including coupled channels. 



The derivations in Section |III| imply that the energy 
corrections should have the same exponential form and 
functional dependence on the radius L at which the wave 
function is zero, independent of the potential and for any 
bound state (although the relationship between L and the 
oscillator determined Li is energy dependent). However, 
there are corrections to Eq. ( [42] ) that become increasingly 
important if L is not sufficiently large. Equation (44) in- 



corporates one such correction but we also have beyond- 



linear energy corrections and the third term in Eq. (27) 



Here we make some representative tests of a direct fit of 
Eq. ([!]) in comparison to applying Eqs. (42) and (44). 

Figure [12] shows results for a quartic potential with a 
moderate depth. The fit to Eq. ([!]) is very good over a 
large range in L 2 for which the energy changes by 30%, 
and the prediction for is accurate to 0.2%. How- 
ever, the fit value of is 1.61 compared to the exact 
value of 1.42. The dashed curve shows the prediction 
from Eq. ( 42 ) using the exact and 700 . It is evident 



that the approximation is very good above Li > 2 but 
increasingly deviates at smaller Li . The modified energy 
correction from Eq. ( 44 ) (dot-dashed curve) matches the 



energy results at the same level as the fit. 



In Fig. 13 examples are shown for square well and 
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FIG. 12: (color online) Energy versus L2 for a quartic po- 
tential well Eq. (|6j for a wide range of N and ftfi (circles) 
(h = H = R = 1). The solid line is a fit to Eq. Q with 
A, koo and E x as fit parameters while the dashed and dot- 

43}. The 
The in- 



dashed lines are predictions from Eqs. \42\ and Eq. 

-1.0115 



horizontal line is the exact energy, E x = 
set illustrates the calculation of the asymptotic normalization 
coefficient (ANC) from the (normalized) wave function. 



FIG. 14: (color online) Energy versus L2 for the deeply bound 
ground state of a Gaussian potential for a wide range of A^ 
and h£l (circles) (h = u = R = 1). These ar e co mpared to the 



predictions of Eq. ( 42 1 (dashed) and Eq. ( 44 1 (dot-dashed) . 
The solid line is a fit to Eq. |l| with A, k^ and Eoo as fit 
parameters. The horizontal dotted line is the exact energy, 
-Boo = -4.2806. 
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FIG. 13: (color online) Energy versus L2 for moderate-depth 
(a) square well Eq. ([3} and for (b) Gaussian potential well 
Eq. ([5} (H = fj, = R = 1) for a wide range of A^ and fif2 (cir- 
cles). The solid line is a fit to Eq. Q with A, k^ and 
as fit parameters while the dashed and dot-dashed lines are 
predictions from Eqs. (42 1 and Eq. ( |44[ ). The horizontal dot- 
ted lines are the exact energies; square well: E^ — —1.5088, 
Gaussian well: E^ = —1.2717 



Gaussian potentials with a moderate depth. Again we 
find a good fit to an exponential fall-off in L 2 , but in 
these cases not only are the energies well predicted (again 
to better than 0.2%) but the fit values of fcoo are within 
5% of the exact results. However, the prediction from 



Eq. (44) actually degrades the agreement for the Gaus- 



sian well compared to the prediction from Eq. (42 ). Fur- 



ther investigation in these cases reveals that the contri- 



keepmg 



42h and (44 1 can fail 



butions from the second and third terms in Eq. (27 1 are 
of comparable size and opposite sign. Therefore, 
only one of them is counterproductive. 

For deeply bound states, Eqs. 
for a different reason. The error in Eq. (42) is propor- 
tional to e~ 4fcooi , so one might expect that the predic- 
tion to become increasingly accurate as the state becomes 
more bound. However, as seen in Figs. [T4| and fl5| results 
for deep Gaussian and exponential potential wells do not 
match this expectation. In deriving the energy correc- 
tions we used the asymptotic form of the wave functions. 
This is valid only in the region r » JJ, where R is the 
range of the potential. The potentials at the smaller val- 
ues of L2 shown in the figures are not negligible. Indeed, 
it is evident from the insets in Figs. [14] and [15] that we 
are not in the asymptotic region for those values of L. 
The lesson is that when applying the IR extrapolation 
schemes discussed in the present paper we need to make 
sure that the two conditions for its applicability are ful- 
filled. First, we need N sufficiently large for L 2 to be the 
correct box size (see Table [l]) . Second we need L 2 to be 
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FIG. 15: (color online) Energy versus L2 for the deeply bound 
ground state of an exponential potential well for a wide range 
of N and HQ (circles) (h = \i = R = 1). These are compared 
to the predictions of Eq. ( [42] ) (dashed) and Eq. p4| (dot- 
dashed). The solid line is a fit to Eq. |l]) with A, koo and 
Eoo as fit parameters. The horizontal dotted line is the exact 
energy, Eoo — —3.3121. 



the largest length scale in the problem under considera- 
tion. 

The results in Ref. [2J and the figures so far are for the 
ground state of the potential. However, the linear energy 
approximation and the specific derivations in the last sec- 
tion should also hold for excited states. This is so because 
the generalization of the results in Subsection |II B| shows 
that (jn / L2) 2 is a very good approximation to the 7 th 
eigenvalue of the operator p 2 for j -C N. In Fig. [l6] 
representative results for excited states from two model 
potentials are shown. We find the same systematics as 
with the ground-state results: the exponential fit works 
very well but the extracted k x is only correct at about 
the 10% level. In assessing the success of Eqs. (42) and 



(44), we note that these excited states in deep potentials 



are comparable to the ground states in moderate-depth 
potentials shown in Fig. 13 The discussion there applies 



here as well, namely that contributions from the second 
and third terms in Eq. (27) are of comparable size and 



opposite sign, so that Eq. (42) alone is a better approxi- 
mation. 

In summary, our tests confirm the expectation from 



Section III that the exponential form of corrections for 
finite HO basis size is universal for different potentials 
and also excited states (and also in one dimension, not 
shown). The leading-order expression Eq. ([I]) is moder- 
ately successful but not quantitative if exact values for 
koo and 7oo are used. This implies that one should not ex- 
pect to accurately extract k^ from a fit to Eq. ([I]). The 



FIG. 16: (color online) Energy versus L2 for the first excited 
states of deep (a) Gaussian Eq. |5| and (b) quartic Eq. Q 
potential wells for a wide range of N and Ml (circles) (h = 
fj, — R = 1). The solid line is a fit to Eq. ([I]) with A, k x and 
Eoo as fit parameters while the dashed and dot-dashed lines 
are predictions from Eqs. (42 1 and Eq. (|44l). The horizontal 



dotted lines are the exact energies for the first excited states; 
Gaussian well: Eoo = —1.2147, quartic well: Eoo = —1.8236 



ment for deep potentials because it is not the dominant 
subleading correction, but we expect it to be the most 
important correction for shallow bound states (including 
the deuteron), which we consider next. 



B. Shallow bound states 

The case of weakly bound states is of special inter- 
est. Figure 17 (a) shows ground-state energies for many 
different N and Ml versus L 2 using Gaussian model po- 
tentials whose parameters are chosen so that the energies 
are the same as the deuteron binding energy (scaled to 
units with H = 1, fj, = 1, R = 1), The prediction Eq. Q 
fails to reproduce the data except at the highest values 
of L2. However, when the correction from Eq. ( |44[ | is 
added there is significant improvement. We also note 
that, contrary to the situation with Figs. (13] and [To] the 
correction from the third term in Eq. ( 27 ) is much smaller 



and of the same sign as the contribution from the second 
term included in Eq. (44). This is consistent with the 



modified energy correction Eq. ( 44 ) is not an improve- 



dot-dashed lines falling below the calculated energies at 
the smallest L 2 values. In Fig. [l7] (b) the same exercise 
is repeated with a model square well. The energies in 
this case are obtained by solving the Schrodinger equa- 
tion exactly with a Dirichlct boundary condition on wave 
functions at r = L. Similar comments as for the model 
Gaussian potential well also apply here. 
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FIG. 17: (color online) (a) Ground-state energy versus Li 
for model Gaussian potential, (b) Energy versus L for the 
square well. The energies for the square well are from solving 
the Schrodinger equation exactly with a Dirichlet boundary 
condition on wave functions at r = L. The dashed and dot- 
dashed lines are predictions from Eqs. ( 42 1 and ( 44 1 . The 



depths of these model potentials are chosen so that the scaled 
energies (with K = fx = R = 1) are the same as the deuteron 
binding energy. 



In Fig. [18] we show analogous results from the deuteron 
calculated with the chiral EFT potential of Ref. [TJ. As 
in Fig. [l7| the modified IR correction Eq. ( 44 ) (evaluated 
using the s-wave ANC) is a significant improvement over 
Eq. (42 1, falling slightly below the calculations at the 



lowest Li values. 



C. Effect of SRG evolution 

As a final test of the universal applicability of the cor- 
rection formulas Eqs. (42 1 and (44 1, we consider a se- 



quence of unitarily equivalent potentials for the deuteron. 
In particular, we use the similarity renormalization group 
(SRG) [26] to evolve the initial Entem-Machleidt po- 
tential to four values of the SRG evolution parameter 
A. Because the transformation is exactly unitary (up 
to very small numerical errors) at the two-body level, 
the measurable quantities such as phase shifts, bound- 
state energies, and ANCs are unchanged. As A decreases, 
the SRG systematically reduces the coupling between 
high-momentum and low-momentum potential matrix el- 
ements, thereby lowering the effective UV cutoff. Thus 
these potentials are useful tools to assess the role of UV 
corrections. 

We first consider results with N and Ml chosen to en- 
sure small UV corrections, as in all prior figures. All the 
quantities on the RHS of formula Eq. ( 44 1 are invariant 
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FIG. 18: (color online) Deuteron energy versus L2 for the 
potential of Ref. [7] . To eliminate the UV contamination we 
only keep points for which Ml > 49. The dashed and dot- 
dashed lines are predictions from Eqs. ( 42 1 and ( 44 1 . The 



horizontal dotted line is the deuteron binding energy. 



under SRG evolution. Therefore, if it is an accurate rep- 
resentation of the IR energy corrections from truncating 
the HO basis, then the EiL-i) vs L2 points for different 
SRG A should lie on the same curve. Figure [19] shows 
that this is the case, and the curve is the same as for the 



unevolved potential in Fig. 18 (Only selected points are 
plotted for readability.) 

Finally, in Fig. [20] we relax the condition that the UV 
corrections are small compared to IR corrections. In par- 
ticular, we fix N at 8 and 12 and scan through the full 
range of Ml. We observe that with increasing L2, each 
of the curves with a given A eventually deviates from 
the universal curve, first with A = 3.0 fnV 1 and then 
later with decreasing A or with higher N . We can un- 
derstand this in terms of the behavior of the induced UV 
cutoff. For fixed N, Eq. tells us that increasing Lo 
means increasing b (or decreasing Ml). But at fixed N, 
A-uv oc 1/6, so the UV cutoff will be decreasing and the 
corresponding UV energy correction increasing. Thus the 
curves at fixed A correspond to the curves seen in conven- 
tional plots of energy versus Ml (e.g., see Ref. [8]). The 
softer potentials (lower A) will have lower intrinsic UV 
cutoffs and therefore they are only affected for larger L 2 . 
The minima for each A are when IR and UV corrections 
are roughly equal. 



V. SUMMARY AND OUTLOOK 

In this paper, we revisited the infrared (IR) correction 
formula derived in Ref. [2] for a truncated harmonic oscil- 
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FIG. 19: (color online) Deuteron energy versus L2 for the 
potential of Ref. [7] evolved by the SRG to four different res- 
olutions (specified by A). To eliminate the UV contamination 
we only keep points for which ftfi > 40. The dashed and 
dot-dashed lines are predictions from Eqs. (42 1 and (441. The 



horizontal dotted line is the deuteron binding energy. 



FIG. 21: (color online) Triton energy versus L2 (here calcu- 
lated with the deuteron-neutron reduced mass) for the two- 
and three-nucleon potential in Ref. [27] unitarily evolved by 
the SRG to four different resolutions (specified by A) with the 
same binding energy |27l 128] . Only larger hfl points are plot- 
ted to minimize the UV contamination. The horizontal dotted 
line is the exact triton binding energy for this interaction. 
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FIG. 20: (color online) The same SRG-evolved potentials as in 
Fig. 19 are used to generate energies, but with N fixed at (a) 
8 and (b) 12 and no restriction on hfl. Thus UV corrections 
are not negligible everywhere. The dashed and dot-dashed 
lines are predictions from Eqs. ( |42[ ) and ( |44[ ). The horizontal 
dotted line is the deuteron binding energy. 



lator (HO) basis expansion, using the simplified case of a 
two-particle system as a controlled theoretical laboratory. 
We used simple model potentials and the deuteron calcu- 
lated with realistic potentials to extend and improve the 
IR formula. We demonstrated analytically that the spec- 
trum of the squared momentum operator p 2 in a finite os- 
cillator basis is identical to the one in a spherical box with 
a hard wall. The minimum eigenvalue of p 2 is (7rfi,/L 2 ) 2 j 
and this identifies L2 as the box radius. While these 
results have been obtained in finite but large oscillator 
spaces, they also hold in practical applications in much 
smaller spaces. We showed how errors parametrized in 
terms of an effective hard- wall radius L from different N 
and HQ, combinations all lie on the same curve, but only if 
the UV error is sufficiently small and, for smaller TV, only 
if L is defined as L2 (see Eq. Q). The determination of 
L2 as the box radius also allows us to extract phase shifts 
from the positive-energy solutions in the oscillator basis. 

The fall-off with L2 of the IR correction to bound- 
state energies is found to be an exponential independent 
of the potential or whether a ground or excited states (or 
whether we are in one or three dimensions). This conclu- 
sion is validated by the derivation and testing of explicit 
formulas for the energy corrections that depend only on 
on measurable bound-state properties: the energy and 
residue of the bound-state pole of the S matrix (or the 
binding momentum and asymptotic normalization con- 
stant). 

Tests on larger nuclei have validated the exponential 
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form Eq. ([I]) with the decay parameter hoc, in the more 
general case associated with the lowest breakup thresh- 
old. Preliminary tests show that L 2 is also the pre- 
ferred definition of L. An example is shown in Fig. |21| 
where triton energies for a two- plus three-nucleon poten- 
tial evolved to four different SRG A (see Refs. [23 [28]) 
lie on the same curve when L 2 is used. A naive fit to 



Eq. ( 42 1 to the triton assuming a break-up into deuteron 
plus neutron yields a binding momentum w 91 MeV 
(jfc^pt = 88.13 MeV) and ANC 7oo rj3 fin" 1 / 2 . The ANC 
is not in agreement with data and previous computations 
where 7oo « 2 fm -1 / 2 was reported [Ml [3D], and sug- 
gests that a more sophisticated analysis is necessary for 
the three-body problem (see also Refs. [3"Tti33] ) . While 
we expect from general considerations that the parame- 
ters of universal curves such as in Fig. [21] are determined 
by asymptotic (and therefore measurable) quantities, it 
remains to be investigated whether simple formulas are 
possible (and whether ANC's might be approximately 
extracted from fits). 

In most of our investigations here we used our ability 
to calculate with very large Ml and N for two-particle 
systems to ensure that the effective UV cutoff was large 
enough to make the UV corrections negligible compared 
to the IR corrections. However, in realistic calculations 
we will not (always) have this luxury. The effects of non- 
negligible UV corrections were shown in Figure [20} By 
working on the other side of the minimum we can isolate 



the UV systematics. Analogous studies to those here but 
on the UV side show that Auv is an appropriate variable 
for the energy correction, but the behavior is not univer- 
sal in the same sense we have identified here. For exam- 
ple, considering different model potentials, ground state 
vs. excited state, and three dimensions vs. one dimension, 
we find there are different functional dependencies (see 
also Ref. [20] )• While some systematic behavior has been 
identified for SRG-evolved potentials 0, further work is 
needed to go beyond the basic form used to make fits. 
Work in this direction is in progress. 
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